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The rigorous application of the magneto-ionic theory to the calculation of reflection 
coefficients for a sharply bounded ionosphere model is carried out. The paper is illustrated 
with computations applicable to the Z)-region or the Z?-region of the ionosphere. The quasi- 
longitudinal approximation is derived from this theory and the range of validity of this 
approximation is illustrated. The restrictions imposed by the use of a sharply bounded 
model ionosphere are discussed. 

1. Introduction 

The classical magneto-ionic theory introduced by Appleton [l], 1 Hartree [2,3], and Booker 
[4] describes the propagation of waves in a homogeneous, anisotropic ionosphere. Various 
methods have been proposed by Bremmer [5], Yabroff [6], and Budden [7] to calculate iono- 
sphere reflection coefficients. The "quasi-longitudinal" (Q-L) approximation [7] provides 
considerable simplification of the calculation. Wait et al. [8] calculated vlf reflection coeffi- 
cients with the aid of this mathematical approximation. 

In this paper, the rigorous theory based on Maxwell's equations for a model ionosphere 
in which the electron density and collision frequency increase abruptly to a finite value at a 
particular altitude 4 and so remain uniformly for all greater altitudes, is reduced to computa- 
tional form. In this paper such a model is called a "sharply-bounded" model and represents 
a "first approximation" of the actual electron density-altitude, collision frequency-altitude 
profiles. The sharply-bounded ionosphere model is a physical approximation at low frequencies ; 
the quality of the approximation decreasing with increasing frequency. The closeness of the 
approximation is also dependent upon the form of the ionosphere electron density-altitude 
profile near the assumed sharp boundary. The theory is illustrated by the calculation of 
certain reflection coefficients applicable to the P-region or the ^-region of the ionosphere. 
Calculations are also presented to illustrate the range of validity of the Q-L approximation. 

2. Theory 

The x?/-plane, figure 1, describes the boundary of the model ionosphere. The vertical 
direction, z describes the normal to the boundary. The region above the #?/-plane (2>0) is 
characterized by a uniform electron density N and collision frequency v. The region below the 
xy-pl&ne (z<C0) is a vacuum. A plane wave is incident on the boundary such that the normal 
to the wave front assumes an angle 4>t (angle of incidence) with the vertical z, and the plane of 
incidence, which can contain either electric E or magnetic H vector, is oriented at a magnetic 
azimuth <p a i.e., the earth's magnetic field vector II m is contained in the yz-pleme, figure 1, with 
a magnetic inclination or dip /. The local coordinate system for the wave is described by 
x', y', z' such that the coordinates x, y, z and x', ?/', z' coincide when the angle of incidence 
<t>i and the magnetic azimuth angle </> a vanish ( <£* = (), <t> a =0). An incident plane wave field E t 
which varies harmonically in time t at a frequency /=co/27r is incident upon the xy-phme, 2 

E t =\Ei\ exp-< i ot — (x sin <fo sin 4>a+y sin <j> t cos <£ a +2 cos <j> t \> • (1) 



1 Figures in brackets indicate the literature references at the end of this paper. 

2 c=speed of light, c~2.998 (108) m/sec. 
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Figure 1. Coordinate systems. 

A resultant wave transmitted into the ionosphere model (z^>0) is then assumed to have the 
form, 



E t = \E t \ exp< i ut — (x sin <j> t sin <j> a +y sin <j> t cos a +2f) V) 



(2) 



where f is, in general, a complex number. Such a wave is described by an electric E and a 
magnetic H intensity vector in Maxwell's classical equations, 3 



VX#+Mo^=0, 



(3) 



VX#-J-eo^=0, 



(4) 



where the permeability and dielectric constant of space are represented by /x and^o, respectively, 
and where the electron convection currents are represented by the vector J—NeV for N 
electrons per cubic centimeter in which each electron has a charge e and travels at a vector 
velocity V; and the equation of motion of an electron, 



m ^+mvV+w(VxHJ+eE=0, 



(5) 



where m is the mass of the electron and t£ m is the earth's magnetic field intensity vector, describe 
the propagation in the ionosphere model. The_ simultaneous solution of these eqs_(3), (4), 
(5) for a wave transmitted into the ionosphere E t (2) upon elimination of the vectors H and V 
can be expressed as the matrix 4 equation, 



A x Ay A z 




~E X 1 


B x By B z 




Ey 


_ ^x ^y ^z _ 




-E 2 _ 


A.= l— at — t* 




s 



=0. 



(6) 



(7) 



3 See for example, H. Bremmer [5] pp. 278 to 280. 

4 See for example, Yabroff [6]. 
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A v =a L a T -i^^, (8) B^a^+^^ry (12) 

A^arf+i-J^, (9) C z =a T $-i-J^, (13) 

B x =a L a r +i-^, (10) C *= a it+ 8 $%) > ( 14 ) 



~s(s 2 -h 2 )' [ } z ~ s(s 2 -h 2 )' 



B y =\-?-a 2 T - , 2 "' , (11) C z =\-d 2 - oU2 "* , (15) 



where, 5 defining the critical frequency squared u cr 2 =Ne 2 lnm and the gy rof requency oo L =/jL H m /m y 

^fi_iZ"l (16) **-*«»/, (17b) 

a£=sin f cos a , (18) 



&=— I' ^ a T = Sln ** SU1 *«> ( 19 ) 

a=sin <£,■. (20) 
h L =—h sin I, (17a) 

Rearrangement of the matrix eq (6) in powers of f results in a quartic equation in f with com- 
plex coefficients, 6 

^r 4 +a 3 r 3 +« 2 f 2 +a 1 f+a =0, (21) 
where, 

«.= («.!« *,_!)* [l--^_] + ( 8 in* a-1) [I+^=| + -^] + _^_, (22) 

*-{ 2 t^^+^&y} (sin ' ^-^+^i)+£|4 < 24 > 



s{s 2 -K 2 ) 



<* 3 =2 -^^-=- ai sec 2 fc, (25) 



a « =1 -, (,>_£!)■ (26) 

3. Determination of the Roots of the Quartic 

The four complex roots of the quartic (21) were found by the Muller [9] iterative method 
applicable to n roots of the nth degree equation, 

Ftf)=a n r+a n -ir- 1 + • • • +a =0, (27) 

where a n , a n -i, • • • a are, in general, complex numbers and a w 5^0. This method does not re- 
quire the calculation of derivatives and necessitates only one value of the function per iteration. 

^ co fr 2=3.18 (10 9 ) N. 

W x,=1.76 (10 7 ) H m , where H m = \H m \ which by convention in geophysical data is called the earth's magnetic intensity in gauss (symbol: r) or 
the earth's magnetic intensity in gammas (symbol: 7) where H m =W y. = Y. 

/c = l/C 2 M0=eo. 

e This equation is equivalent to equations given by Bremmer [5] eq (18), p. 291 and YabrofT [6] p. 751. 
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Each iteration in the process of rinding one root is obtained from a calculation of the nearer 
root of the quadratic equation which passes through the last three points of the function, 
F(£). This quadratic equation is, in general, complex, i.e., it possesses complex coefficients and 
complex roots. A LaGrange interpolation formula was employed in the process of finding the 
root. The function, /(f), in the LaGrange interpolation formulas, 

M)=A n r+A n _ 1 r~ 1 + . . • +A , (28) 

is an nth. degree equation with the same ordinates as F({) for n + 1 distinct abscissas [10] as 
f», tn-i • • • fo- The LaGrange interpolation formula through the three points; f*_ 2 , F($ t - 2 )) 

U-u F{U-i)) and f f , F{U) of F({) is, 

J(t)=M 2 +A^+A 0) (29) 

where A 2 , A u and A are evaluated from the requirement: 

Upon solving for A 2 , A u A and using the quantities 0=f— £«, Pi=£i— f*-i, 0«-i=?«-i— £ f _ 2 , 
\=P/p i) Xi=/3j//3 z _i, and 5^=1 +X<, the LaGrange interpolation formula, eq (29), becomes a 
quadratic in X: 

+\d- 1 [F(^ 2 )\ 2 i -F(^ 1 )8]+F(^ i ) fa+BM+Ffa). (30) 
f f+1 is found for the condition, /(f) =0, solving for X and employing the relationship, 

X=X W =£±^. (31) 

hi Si-1 



Rationalizing the numerator of the standard quadratic formula, and solving for X z+1 , 



(32) 



.! -<*> 



gfjgi-mti) fat [*Xr<-2) Xz-^(fz--i) «i+/^(fi)] 

where, 

^=^(f,_ 2 )X?-^(f < _ 1 )^+^(f,)(X i +5,). (33) 

The sign choice in the denominator of \ i+1 (32) is resolved by selecting the value for the larger 
denominator. This choice of \ i+l gives £ <+1 the root of the LaGrange interpolation formula 
(29) which is nearer f*. The initial values of f and F(£) can be taken as follows: f = — 1, 
i^(f )=a 2 — ai+a , fi=l, -F(fi)=a 2 +a 1 +a , f2 = 0, and F(^ 2 )=a . The iterative process is 
terminated with f*, the derived root when, 

If* Tz-i 

IM 

where e is a predetermined number. 

Upon determination of a root, the degree of the polynomial can be reduced by dividing 
F(£) by the root. The new coefficients can be found by, 

a{=f w a{. 1 +a < ,i=0,l,2 . . ., 

with a_i = 0, where a\ replaces a% and f w is the root already calculated. 

Two pairs of roots, f , where each root represents either an "ordinary" or an "extraordinary' ' 
wave propagated in the ionosphere, can be identified as upgoing {-\-z direction, fig. 1, usually a 
fourth quadrant 7 f) and downgoing ( — z direction, fig. 1, usually a second quadrant f) waves. 
The upgoing pair determine the reflection coefficient for waves incident on the ionosphere from 
the region below the ay-plane (2<0). 

7 It is also possible to find roots in the 3d quadrant under certain conditions. See for example, Booker [4]. 
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The quartic (21) is satisfied by the classical Appleton-Hartree formula, 8 

2 (i-i) 

7,2=1 ; > 4 s/ => (34) 

2 ( i -4)-? sin2 ^ w sin4 * + ^ ( 1- *) 2 <:os2 * 

where, 

fW-sin 2 ^, (35) 

and, 

f 1 

cos ^= — sill 7H — sin 6< cos fl cos /, (36) 

v v 

sin ^=- [rj 2 — f 2 sin 2 /—sin </>* cos </> (sin <£*• cos </> a cos 2 1— 2f sin / cos 7)] 1/2 . (37) 

Although the Appleton-Hartree formula (34) provides a solution of the quartic for the up- 
going waves (the dual sign in the denominator (34) represents ordinary and extraordinary 
waves), it is also necessarv to satisfy the auxiliary eqs (35) to (37). But the iterative solution 
(27) to (33) is identical with such a process (34) to (37). 

4. Determination of the Reflection Coefficients 

The reflection coefficients are determined by five boundary conditions which express the 
principle of continuity of the tangential E and //fields and the normal // field 9 at the boundary 
(a^/-plane, fig. 1) of the model ionosphere. The E x field, figure 1, immediately above and im- 
mediately below the boundary can be equated, 

E x , t cos <t> a +E yfi cos <j>i sin <j> a -\-E x , T cos <\> a —E y , T cos 4> t sin (f)a=QoE yo +Q e E ye ) (38a) 

also the E y field, 

—E x ,i sin 4> a +E yfi cos </>*• cos <f> a —E x , f sin $ a —E y , T cos cj> t cos cj) a =E y0 +E ye ; (38b) 

the H x field, 

E x , t cos </>* sin 4> a —E yfi cos <t> a —E x , r cos 4> t sin cf) a —E y , r cos cj) a = (a L P —{ )E y0 + (a L P e —^ e )E ye ) 

(38c) 
the H y field, 

E Z f t cos <j>i cos 4> a -\-E y ,i sin <\> a —E x t r cos 0* cos <t> a -\-E v r T sin <fr a 

= ttoQo-a T Po)E yo +tt e Q e -a T P e )E ye ; (38d) 
and the H z field, 

—E xn sin <t> t —E x , r sin (j> i ={aT—a L Q )E yo J r{a T —a L Q e )E ye ) (38c) 

where 10 (see fig. 1) i?^ is the electric field vector normal to the plane of incidence of the incident 
wave, E v n is the electric field vector in the plane of incidence of the incident wave, E x , r is the 
electric field vector normal to the plane of incidence for the reflected wave, E y , T is the electric 
field vector in the plane of incidence for the reflected wave, also, 



8 See for example, H. Bremmer [5], p. 282, eq (8). The direction \f/ is complex and double-valued (\f/ , ^e) as a result of the two roots, f and 
f«. The authors do not concur with a certain previous interpretation of the Appleton-Hartree formula, op. cit. [6], p. 751. It can be shown that 
the classical Appleton-Hartree formula, eqs (34) to (37), after considerable ado reduces precisely to the quartic, eq (21), the coefficients of which 
have been described, eqs (22) to (26); the latter of which are by definition valid for both zero 0>=0) and finite (»*>0) collision frequency. The 
authors must therefore conclude that the classical Appleton-Hartree formula is valid for finite collision frequency (v>0). 

9 See for example, Stratton [11], p. 483. The fifth boundary condition (38e) expresses the continuity of the normal H field. As a consequence 
of this condition (38e) it can be shown that (see eqs (41a-h)), 

-T mm -[(a T -a L Qo) U mo +(a T -aL Q.) I7«j/sin «/»,— 1=0. 

The fifth condition (38e) must be satisfied automatically by the previous four conditions (38a-d), and hence the above expression (derived from 
38e) can be employed as an independent partial check on the entire computation. 

10 The subscripts o and e refer somewhat arbitrarily to the two roots off corresponding to "ordinary" and "extraordinary" waves. 
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n— Ex 



(39) 



Ik 

Ey 



(40) 



Four reflection coefficients and four " transmission' ' coefficients can now be defined, 11 



■*- ee T7i y 



T =■ 

-*- em, 



T = 



E x t r 

Eyfi 

E y r T 

_ E x t T 
Ez't 



(41a) 



(41b) 



TT Ey 

U eo J7i f 



u m 



_Eyo_ 

E x , t 



Ey, 



(41c) U ee =^-, 



(41d) 



77 ye 

*-J me — 77i f 



(41e) 



(41f) 



(41g) 



(41 h) 



where, T ee refers to the plane wave ionosphere reflection coefficients with electric, E, vector in 
the plane of incidence of the incident wave and the electric, E, vector in the plane of incidence 
for the reflected wave; T em refers to the same incident wave but the magnetic component 
(abnormal component) is in the plane of incidence for the reflected wave. Similar normal 
(electric) and abnormal (magnetic) components arise in the reflected wave resulting from an 
incident wave with magnetic component in the plane of incidence which components are de- 
scribed by the coefficients, T me and T mm , respectively. The quantities, U eo , U mo , U ee , and U me 
define the nature of an ordinary wave (second subscript o) and an extraordinary wave (second 
subscript e) excited in the model ionosphere by an incident wave with the electric vector in 
the plane of incidence (first subscript e) or magnetic vector in the plane of incidence (first 
subscript m) of the incident wave. 

The matrix solution of the boundary conditions (38a, b, c, d) for an incident wave with the 
electric, E, vector in plane of incidence, 



a, 


a 2 


^3 


a 4 




T 

-*- ee 




^oe 




b 2 

c 2 




b 4 

c 4 




T 

-l em 

u eo 


+ 


b e 
^oe 


dr 


d 2 


(h 


di_ 




.u ee _ 




_ Moe _ 



=0, 



(42) 



and the matrix solution for an incident wave with the magnetic, H, vector in the plane of 
incidence, 



(43) 



completely define the four reflection coefficients: T ee , T em , T mm , T me ; the " transmission" co- 
efficients: U eo , U ee , U m0} U me , although of secondary interest, are a byproduct of the solution, 
where, 



'&! 


a 2 


a z 


aC 




FT 1 

- 1 - me 




^om 


6i 

Ci 


c 2 








T 

-*- mm 
U mo 


+ 


V om 


d 1 


d 2 


d, 


6? 4 




U me 




(lorn 



a oe =cos fa sin a , 
a m=zos 4> a , 
b oe =cos <j> t cos </> a} 



(44a) 


b om =— sin <f> a , 


(44b) 


C oe =—COS (j) a , 


(44c) 


c om =cos <j> t sin 4> a 



(44d) 
(44e) 
(44f) 



11 Bremmer's textbook notation for the reflection coefficient, T, is followed in this paper [5] pp. 286 to 295. Budden [7] and Wait [8] have em- 
ployed for the reflection coefficients (41a, b, c, d) the notation \\R\\, \\Ri, iR\\, jlRi, respectively. However, for application to multihop geo- 
metric-optics (Bremmer [5] p. 194 et seq.) the " T" notation can be more readily discerned in the complicated formulas containing both ionosphere 
and ground reflection coefficients, in which the symbol R, has already been used for the latter. 
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(44g) 


b±=-l, 


(44p) 


(44h) 


Ci = — cos </>i sin fl , 


(44q) 


(44i) 


c 2 =— COS </>„, 


(44r) 


(44j) 


^3 = fo #L* oj 


(44s) 


(44k) 


c 4~ f e &Lle> 


(44t) 


(441) 


di = — cos 2 - cos a , 


(44u) 


(44m) 


d 2 =sin a , 


(44v) 


(44n) 


d 3 =a T P — toQo, 


(44\v) 


(44o) 


(14 = (It1 e £eQe' 


(44x) 



d oe =sin <j) a , 

d om = COS 4>i COS (f> a , 

ai=cos a , 

a 2 = — cos Z - sin <£ a , 

#4~ Vej 

bi=— sm<t> a) 

b 2 = — cos 0< cos c/> a , 

63= -1, 

The solution of the matrices (42 and 43) was obtained by an application of Crout's [12] 
modified Gaussian method. Although it is quite possible to solve the matrix analytically and 
there are several numerical methods for the solution of matrices, the Crout method was em- 
ployed because of its adaptability to electronic data processing and the requirement that the 
method chosen must be readily applicable to the solution of matrices with complex coefficients 
of the variables. The Crout method is applicable to n linear equations in n unknowns with 
either real or complex variables. The matrices (42) and (43) are of the form, 

&ll#l T #12#2 T • • • ~\&\n%n = c l 
#21«^1 ~T~ #22«^2 I ■ • • \Q / 2n$n == C2 



#»l2a~"r#n2#2~T • • • \&nn — C m 



(45) 



where the n 2 coefficients a tj and the n right-hand numbers c t are given (44a-x). Since it is neces- 
sary to know the position of each element in the matrix to develop the method of solution, 
each element has a double subscript, that is, element a tj is the element in the ith row and jth 
column. 

The augmented or given matrix, ||a r7 || for the system (45) is formed by adjoining the 
right-hand numbers (the c's) to the matrix formed by the coefficients of the unknowns. This 
matrix can be expressed in rectangular array, 







#110-12 ' 


. . a ln 


Ci 


a 


u\\ 


= #21#22 • 


. . (l2n. 


C2- 






^nl^nl • 


• ("Tin 


Cji 



(46) 



First the elements of the auxiliary matrix, ||a^|| are determined, 



dll#12 • • • din 



\&ij\\ — #21#22 • • • #2« 



(47) 



&wl#«2 



&nn \Cn 



From this matrix a final column matrix, \\Xi\\, with the elements for the required variables 
Xi . . . , x n can be obtained, 
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\Xi\\ — #2« 



(48) 



Each element in both the auxiliary and final matrix is evaluated from previously calculated 
data by a sequential operation. 

The following equations were used to calculate the elements a is from the elements a a of 
the given matrix: 

i-i 
a>ij=a>ij—H aMp (i>j) (49) 



da L k-i J 

,_1T t\ , . 1 

Ci — ~T Ci 2—1 a ik a kj V 

&u L fc=i J 



(50) 



(50a) 



where any summation whose lower limit exceeds the upper is zero. Since each a[j in the auxiliary 
matrix, except those in the first row and first column, depends upon the value of previously 
calculated elements in the same matrix, the elements must be found in a certain order. The 
elements of the first column are found (49) as follows: 

a n =(Lu, a 2 i = a 2 \, a 3 i=a f su and a^=a 4 i. 

The remaining elements of the first row are (50), 



a 12 =^r> d\z=---ri du=~7^} and c[- 



012 

a n 



On 



-A. 
~a' n 



With the aid of these values, the remaining elements can be found (49) or (50). a^ can be made 
available for the solution of the matrix ||a^|| by the following sequential calculation: The n 
elements of the first column, the remaining n of the n-\-\ elements of the first row, the re- 
maining ?? — l elements of the second column and second row, the remaining n—2 elements 
of the third column and third row, continuing in this sequence until the elements are calcu- 
lated [9]. 

The values of x are found from x n to x x from the elements of the auxiliary matrix, 



Xi — Ci 2—1 UikXk- 
k=i+l 



(51) 



It is evident (51) that x n =Cn. The remaining solutions x n -i, x n _ 2 . . ■ , x 2 , x x can be found (51) 
with the aid of x n . The computations presented in this paper were checked by substitution of 
Xi, x 2 , . . . x n in each of the original equations. An inductive-type proof of the Crout method 
has also been given by Crout [12] and Hildebrand [10]. 

The analytical expressions for the complex numbers P , P € , Q , Q e can be derived from the 
definitions (39) and (40) with the aid of the matrix (6) resulting from the simultaneous solution 
of Maxwell's eqs (3), (4), and the equation of motion of an electron (5) with the following 
result : 



P= 



Q- 



-b^%] 


[i < f 2 ,._*.] i[ a *°* v_V] 


l a r? V-A 2 J 


r s>-hi - 

L s(s 2 -h 2 )] 


\} < f 2 ,_d [^f+v-A\] 





-[ 



1-a 2 



s 2 -hl 
's(s 2 -h 2 ) 



[ aL a T _i _^_J + ^ rf+i _^L_J 



b-*-M>I 



l-a£-f 2 



(f-ti 



]-[a 3 



r+i 



s 2 -h 2 



ad- 



h L h T 



s(s 2 -h 2 )J 



Orf— i 



li-i 



s 2 -h 2 



(52) 



(53) 
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where the subscripts o and e for P , P e , Q , Q ey refer to the two roots of f which represent 
upgoing waves in the model ionosphere. 

5. The Q-L Approximation 

Budden [7] and Wait [8] have evaluated reflection coefficients for the sharply bounded 
model ionosphere employing the mathematical Q-L approximation. This approximation has 
been employed to interpret experimental measurements of skywaves propagated at vlf between 
the surface of the earth and the ionosphere in the northern hemisphere. 

Such an approximation can be derived from the rigorous formulation presented in this 
paper. The model ionosphere reflection coefficients can be written analytically from the rigor- 
ous matrix eqs (42) and (43), 

T B l A 2 —A l B 2 A .-xVdi hi .-if^ c 2 l ( v 

T "=AsB 2 -B z A 2 > (54) ^ 3=Al U""d" A2 larcj (60) 

B.A-ABs 5 _ A MK^y„ A -iK 41 (61) 

A 2 K Z — J3 2 A 3 \_a A 4 J L a 4 (IaJ 



T = 

-*- em 

T = 

-*- me. 



A3B2 B S A 2 



<« *-*- , &-£h" , [jS-& (62) 

^l I= A r ir^£_^"|_A 2 -T— ~— T (58) d=Ar l [— -M-A.' 1 !"— -— "I 

L«4 t> 4 J L«4 C 4 J L«< &jj L»4 c 4 J 

^= Ari [s4il-^ i [i-3' (59) *M*-fcM*-fc) 



(64) 
(65) 



where, 



Ai-5*!, (66) A 2 =^-A (67) A 3 =^-^ (68) 

a 4 64 a 4 c 4 a 4 r/ 4 

The Q-L approximations are as follows: 

$ a =0, (69) 

Qo 



Qe = 



COS O 

cos 6 e 
or, 

Vo 



(70) 


P =- 


sin </>« 
77 cos 0„ 


(71) 


P e =- 


sin 4> t 
Ve cos 0/ 


(74) 


Q.f. 


= h 



(72) 
(73) 

(75) 



which ((74) and (75)) is approximately true for vlf skywaves reflected in northern hemisphere. 
After considerable ado, the Q-L reflection coefficients with the aid of the approximations (69) 
to (75) and SnelPs law, 12 

sin (j> i =rj 0te sin Q 0>e , (76) 



12 The validity of Snell's law for complex directions, 9 and 0«,i s also implied by eqs (36) and (37). 
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reduce as follows: 

T ee =Ab l {(Vo+Ve) (COS 2 4>i — COS d COS 6 e ) + (rj 7] e — 1) (cos o +cos 0«) COS </>*}, (77) 

r m =A 6 _1 {2 i cos <£, (r) cos O — r? e cos e )}, (78) 

T me =A b - 1 {2 i cos 0« (r) cos e — ?? e cos o )h (79) 

T mm = A^{ (vo+Ve) (cos 2 (f)i— cos O cos e ) — (voVe— 1) (cos O + COS e ) COS <£,}, (80) 

where, 

A 6 =(7? 4-?7 e ) (cos 2 <^+cos O COS0 e ) + O7</>?e+l) (cos O +COS e ) COS CJ) i} (81) 

which formulas agree with those independently derived by Budden [7], The indexes of 
refraction of the model ionosphere, rj and rj e , corresponding to upgoing "ordinary" and "ex- 
traordinary" waves can be calculated rigorously as the roots of the quartic in f (21). However, 
as a result of the matrix asymmetry introduced by the Q-L approximation (74) and (75) an 
interchange of the roots of the quartic (77 and rj e ) results in an ambiguity, ±7r, in the argument 
or phase of the abnormal components, arg T em and arg T me (78) and (79). At very low fre- 
quencies, the indexes of refraction rj and 77 e are chosen such that arg T em and arg T me are close 
to 7T. The values 77 can also be calculated from the Appleton-Hartree formula (34). Budden 
[7] and Wait [8] employ additional approximations to this formula, 13 

V 2 o.e=l~i~exp[±i<l>i], (82) 

in which the plus ( + ) sign before cj> x refers to rj and the minus ( — ) sign refers to 77 e and 
where, 

tan^=^, (83) £— _^L~ (84) 

» « coV" 2 +«| 

6. Computations and Discussion 

The results of a calculation 14 of model ionosphere reflection coefficients which can, under 
certain conditions, represent the reflection of waves from the Z?-region or the jE'-region of the 
ionosphere are illustrated, figures 3 to 32. These data have been presented as a function 
of frequency, for various angles of incidence on the ionosphere, (/>*. The application of these 
reflection coefficients to the geometric-optical theory is quite simple, 15 however, the relation 
between angle of incidence, 4> u distance of observer on the surface of the earth from the source 
on the surface of the earth, d/j (j^order of hop or time-mode, i.e., j=l, 2, 3 . . .), and the 
altitude of the boundary of the model ionosphere have been presented for the convenience 
of the reader, figure 2. The results of a rigorous calculation of the reflection coefficients of 
the model ionosphere are compared with the results of the Q-L approximation, figures 29 to 32. 
Thus, for magnetic inclinations near 60°, the Q-L approximation is quite satisfactory at fre- 
quencies less than 8 kc for the normal components, T ee and T mm . The most serious discrepancy 
occurs in the prediction of the abnormal components, \T em \ (fig. 29, for example). At higher 
vlf frequencies (> 8 kc) and at If, the Q-L theory fails except in certain special cases (|T me |. 
fig. 32, for example). 

The effects of magnetic azimuth </> a , figures 19, 20 (compare with figs. 3, 4) magnetic 
inclination or dip /, figures 21 to 24, magnetic intensity H m , figures 25 to 28, and distance or 
angle of incidence (/>*, figures 3 to 18, on the reflection coefficients are illustrated. It is of interest 
to note that an increase in the earth's magnetic intensity H m increases the abnormal components 
T em , T me , figures 27, 28. It is also of interest to note that the reciprocity theorem does not 

" The validity of these approximations has been discussed by Budden [7]. 

14 The calculations were performed on the IBM-650 electronic data processing machine. 

15 See for example, H. Bremmer [5] p. 153 et seq. 
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Figure 2. Geometric-optical relation between the 
angle of incidence, 4>\, distance from the source, d/j, 
(j = order of skywave hop, ] = 1, 2, 3 . . .) and 
altitude of boundary of model ionosphere, h. 
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Figure 3. Model ionosphere reflection coefficients; 
<f>i = 75.08°, sin 4^ = 0.966856, N = 870, v = 4(10«), 
R m =0.5, 4> & = 0, 1 = 60° (for example, fig. 2, 
d/j = 329 miles, h = 65 k m ) . 
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Figure 4. Model ionosphere reflection coefficients; 
4> { = 75.08°, sin 4>, = 0.966 356, N = 870, i> = 4(10«), 
U m = 0.5, <t>* = 0, 1 = 60° (for example, fig. 2, 
d/j = 329 miles, h = 65 km) . 




10" 10° 

FREQUENCY,f,CpS 

Figure 5. Model ionosphere reflection coefficients; 
4>i = 80.89°, sin 4> { = 0.985948, N = 870, v = 4(10*), 
H m = 0.5, <f> & = 0, 1 = 60° (for example, fig. 2, 
d/j = 621 miles, h = 65 km). 
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Figure 6. Model ionosphere reflection coefficients; 
4>i = 80.89°, sin 4> x = 0.9859 48, N = 870, p = 4(10 Q ), 
H m = 0.5, 4> & = 0, 1 = 60° (for example, fig. 2, 
d/ i = 621 miles, h = 65 km). 
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Figure 7. Model ionosphere reflection coefficients; 
4>i = 81.79°, sin 4>i=0.989753, N = 870, v=4(10*), 
K m =0.5, 4>*=0, 1 = 60° (for example, fig. 2, 
d/j = 1000 miles, h = 65 km). 
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Figure 8. Model ionosphere reflection coefficients; 
<f>i = 81.79°, sin S = 0.9897 '53, N = 870, v = 4(10 6 ), 
U m = 0.5, <j> & = 0, 1 = 60° (for example, fig. 2, 
d/j = 1000 miles, h=65 km). 
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Figure 9. Model ionosphere reflection coefficients; 
4>i = 48.25°, sin 0, = 0.685204, N = 1200, v=10\ 
B. m =0.5, <f> & =0, 1 = 60° (for example, fig. 2, 
&j] = 100 miles, h = 85 km). 
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Figure 10. Model ionosphere reflection coefficients; 
<J>i = 43-25°, sin <f> { = 0.685204, N = 1200, v = 10\ 
H m = 0.5, <f> 2 =0, 1 = 60° (for example, fig. 2, 
dlj = 100 miles, h = 85 km). 
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Figure 11. Model ionosphere reflection coefficients; 
4>i = 71J3°, sin <fc = 0.946258, ^ = 1200, v = 10*, 
H m = 0.5, 4> & = 0, 1 = 60° (for example, fig. 2, 
dj] = 329 miles, h = 85 km). 
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Figure 12. Model ionosphere reflection coefficients; 
<t> [ = 71.13°, sin <f> x = 0.946258, N = 1200, v = l&, 
H m = 0.5, <f> & = 0, 1 = 60° (for example, fig. 2, 
d/j = 329 miles, h = 85 km). 
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Figure 13. Model ionosphere reflection coefficients; 
cf>i = 78.17°, sin ^ = 0.978751, N = 1200, v = 10«, 
B. m = 0.5, <f> a = 0, 1 = 60° (for example, fig. 2, 
djj = 621 miles, h = 85 km). 
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Figure 14. Model ionosphere reflection/ coefficients; 
H<f>i=78. 17°, sin $-, = 0.978751, N = 1200, v = 10\ 

H m = 0.5, $ a = 0, 1 = 60° (for example, fig. 2, 

d/j = 621 miles, h = 85 km). 
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Figure 15. Model ionosphere reflection coefficients; 
$, = 79.64°, sin $-, = 0.983688, N = 1200, v=10\ 
H m =0.5, $ & = 0, 1 = 60° (for example, fig. 2, 
d/j= 800 miles, h = 85 km). 
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Figure 16. Model ionosphere reflection coefficients ; 
$-, = 79.64°, sin $, = 0.983688, N = 1200, v = 10\ 
H m =0.5, $ & = 0, 1 = 60° (for example, fig. 2, 
d/i. = 800 miles, h=85 km). 
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Figube 17. Model ionosphere reflection coefficients; 
$, = 8040°, sin $1 = 0.985987, N = 1200, v = 10\ 
K m =0.6, $ a = 0, 1 = 60° (for example, fig. 2, 
d/j = 1000 miles, h = 85 km). 



1.0 
E* 0-6 





1 ' '■'! ' ! j I ""1 ' . 


I , i '1 ll| 1 1 Mllll 




/ M mml 






/ arg T me ^^^^^ 


- 


or Q T mm 


t ^ 




l T meK^ -*■ ^ 


■ . .-.-.I ........ 



10' 10' 10" 10° 10° 10' 

FREQUENCY.f.cps 

Figure 18. Model ionosphere reflection coefficients ; 
$-, = 8040°, sin $1 = 0.985987, N = 1200, v=10\ 
H m = 0.5, $ a = 0, 1 = 60° (for example, fig. 2, 
d/j = WOO miles, h = 85 km). 
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Figure 19. Model ionosphere reflection coefficients t 
illustrating the effect of magnetic azimuth, $ SL 
(compare with fig. 3); $, = 75.08°, sin $ { = 0.966356 T 
N = 870, v = 4(10«), R m = 0.5, $ & =180°, 1 = 60* 
(for example, fig. 2, d/] = 329 miles, h = 65 km). 
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Figure 20. Model ionosphere reflection coefficients, 
illustrating the effect of magnetic azimuth, <f> & 
(compare with fig. 4); <j> i = 75.08° i sin <f> { = 0.966356, 
N = 870, v = 4(W 6 ), H m = 0.5, <t> & = 180°, 1 = 60° 
(for example, fig. 2, d/j = 329 miles, h=65 km). 
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Figure 21. Model ionosphere reflection coefficients, 

illustrating the effect of magnetic inclination or dip, 

I (compare with figs. 3, 23); 4> l =75.08°, sin <t>i = 

0.966356, N = 870, v = 4(10*), U m = 0.5, <f> & =0, 

1 = 0. 
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Figure 22. Model ionosphere reflection coefficients, 
illustrating the effect of magnetic inclination or dip, 
I (compare with figs. 4, ^4) ; <f>i = 75.08°, sin 
<j>i = 0.966356, N = 870, v = 4(10*), U m = 0.5, 
<t> & =0, 1 = 0. 
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Figure 23. Model ionosphere reflection coefficients, 
illustrating the effect of magnetic inclination or dip, 
I (compare with figs. 3, 21); <j> x =75.08°, sin 
^ = 0.966356, N = 870, v = 4(10 Q ), R m = 0.5, 
<f> a = 0, 1 = 90°. 
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Figure 24. Model ionosphere reflection coefficients, 
illustrating the effect of magnetic inclination or dip, 
I (compare with figs. 4, 22); (f>i = 75.08°, sin 
<^ x = 0.966356, K=870, v = 4(10*), B m = 0.5, 
<f> & = 0, 1 = 90°. . ^ 
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Figure 25. Model ionosphere reflection coefficients, 
illustrating the effect of the intensity, H m , of the 
earth's magnetic field vector, H m (compare with figs. 
3, 27); 4>i = 75.08°, sin cf> { =0. 966356, 7$ = 870, 
v=4(W), R m =0.2, <f> & = 0, 1 = 60°. 



282 




I(T I0 J 

FREQUENCY, f.cps 

Figure 26. Model ionosphere reflection coefficients, 
illustrating the effect of the intensity, H m , of the 
earth's magnetic field vector, H m (compare with 
figs. 4, 28); <f> x =75.08°, sin <f>\ = 0.966356, N = 870, 
v = 4(10*), 11^ = 0.2, <f> & = 0, 1 = 60°. 
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Figure 27. Model ionosphere reflection coefficients, 
illustrating the effect of the intensity, II m , of the 
earth's magnetic field vector, H m (compare with 
figs. 3, 24); d>i = 75.08°, sin { = 0.966356, N = 870, 
v = 4(10*), ll m =l, <j> & = 0, 1 = 60°. 
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Figure 28. Model ionosphere reflection coefficients, 
illustrating the effect of the intensity, H m , of the 
earth's magnetic field vector, H m (compare with 
figs. 4, 25) ^ = 75.08°, sin <f> { = 0.966356, N = 870, 
v=4(W), H m =l, 4> & = 0, 1 = 60°. 
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Figure 29. Comparison of reflection coefficients 
calculated by the rigorous method for the model 
ionosphere with the Q-L approximation; (/>; = 50.84°, 
sin 0j = 0.7? '5353, N = 870, v = 4(10*), R m = 0.5, 
<l>a = 0, 1 = 60° (<f> & and I do not apply to Q-L 
method; also note: fig. 2, d/j = 100 miles, h = 65 
km). 
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Figure 30. Comparison of reflection coefficients 
calculated by the rigorous method for the model 
ionosphere with the Q-L approximation ; 4> x = 50.84°, 
sin 4>\ = 0.775353, ^ = 870, v = 4(10«), ll m = 0.5, 
(j> & = 0, 1 = 60° (<f> & and I do not apply to the Q-L 
method; also note: fig. 2, d/j = 100 miles, h = 65 km) . 
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Figure 31. Comparison of reflection coefficients 
calculated by the rigorous method for the model 
ionosphere with the Q-L approximation; <f>i = 81.37°, 
sin 4,1 = 0.988681, N = 870, v = 4(10*), Ji m = 0.5, 
a = 0, 1 = 60° (0 a and I do not apply to the Q-L 
method; also note: fig. 2, d/] = 800 miles, h = 65 
km). 
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Figure 32. Comparison of reflection coefficients 
calculated by the rigorous method for the model 
ionosphere with the Q-L approximation; <f>\ = 81.37°, 
sin $-, = 0.988681, N = 870, v = 4(10 & ), R m = O.S, 
(j> & = 0, 1 = 60° (cf> a and I do not apply to the Q-L 
method; also note: fig. 2, d/j = 800 miles, h = 65 
km). 



hold, i.e., south-north propagation, <£ a =0, figures 3, 4, is not precisely equivalent to north-south 
propagation, </> a =180°, figures 19, 20. It is important to note that the magnetic dip /, figures 
21 to 24, and_the magnetic azimuth <j> aj figures 19 to 20, i.e., the direction of the earth's magnetic 
field vector H m is a very significant consideration, not only at If but also at vlf. 

Since the actual ionosphere electron 16 density versus altitude profile is quite variable, 
considerable caution is necessary in the application of this model to the interpretation of 
experimental data. It can, in general, be stated that the model is valid when the boundary 
of the ionosphere is sufficiently sharp (relative to a wavelength) at the frequency of the incident 
wave under consideration. The precision with which this model represents the actual 
ionosphere can be considered a first approximation at low frequencies, the degree of approxima- 
tion decreasing with increasing frequency. 

7. Conclusions 

The reflection coefficients for the sharply bounded model ionosphere can be readily 
evaluated with the aid of electronic data processing techniques for the solution of the quartic 
equation and the tensor matrix, thus eliminating the necessity for the Q-L approximation. 
The Q-L theory is a crude mathematical approximation taking account of the intensity bu t neglect- 
ing completely the direction (magnetic azimuth and dip) of the earth's magnetic field vector and 
if used must be applied with considerable caution not only at If but also at vlf frequencies- 
The complete and rigorous magneto-ionic solution for the assumed model ionosphere not only 
eliminates the necessity for such an approximation but also illustrates the application of 
modern numerical analysis techniques which can in many cases obviate such mathematical 
approximations, leaving only the physical approximations implied by the assumed model. 

The model ionosphere presented in this paper represents a first approximation of the 
electron density-altitude profile, the degree of approximation decreasing with increasing 
frequency and the validity of the model dependent upon the sharpness of the actual ionosphere 
boundary relative to a wavelength. The results of this paper suggest a refinement of the 
theory which would take into account a "less sharply bounded' ' model ionosphere. 



[13]. 



i« Values of electron density, N, and collision frequency, v, employed in this paper were estimated from recent data, see for example, Waynick 
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